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Abstract. Fully dynamic simulations of heavy ion collisions at RHIC and 
at LHC energies within the perturbative QCD-based partonic transport model 
BAMPS (Boltzmann Approach to Multi-Parton Scatterings) are presented, 
r^^ ' focusing on the simultaneous investigation of jet quenching and elliptic flow. The 

^1^' model features inelastic 2 3 processes based on the Gunion-Bertsch matrix 

element and has recently been extended to include light quark degrees of freedom, 
allowing for direct comparison to hadronic data on the nuclear modification factor 
via a fragmentation scheme for high-p^ partons. The nuclear modification factor 
of neutral pions in central Au-|-Au collisions at RHIC energy is compared to 
experimental data. Furthermore first results on the nuclear modification factor 
and the integrated elliptic flow of charged hadrons in Pb -|- Pb collisions at LHC 
are presented and compared to recent ALICE data. These investigations are 
complemented by a study on the suppression of D-mesons at LHC based on elastic 
interactions with the medium. 



[-~ ' 1. Introduction 

o . 

The suppression of higli-pr spectra witli respect to a scaled p + p reference, quantified 
in terms of the nuclear modification factor Raa, and the strong collectivity of the 
medium, quantified in terms of the Fourier coefficient V2, the elliptic flow, have 
^ . been established by experiments at the Relativistic Heavy Ion Collider (RHIC). First 

^ ' results from the recently commissioned Large Hadron Collider (LHC) have confirmed 

these findings at an order of magnitude higher collision energy, showing remarkable 
quantitative similarities beyond a good qualitative agreement. The suppression of 
particles with high transverse momentum is commonly attributed to an energy loss 
on the partonic level during the quark-gluon plasma (QGP) stage of the evolution 
of the medium. Comparison of V2 measurements to hydrodynamic simulations have 
established that the elliptic flow builds up early and is thus also governed by the 
evolution of the QGP medium. It is notoriously difficult to combine both aspects of the 
QPG evolution — the suppression of rare, high-pr probes and the collective behavior of 
the bulk particles — into a common model. The work presented in this article explores 
the capabilities of microscopic partonic transport models in this respect. 
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2. The transport model BAMPS 

The microscopic transport model BAMPS (Boltzmann Approach to Multi-Parton 
Scatterings) [1] is applied to simulate the time evolution of the hot partonic medium 
that is created in heavy ion collisions at RHIC and LHC. It is based on matrix elements 
in leading order perturbative QCD (pQCD) and consistently features inelastic 2 4-> 3 
interactions. Rapid thermalization of the simulated partonic medium on the order 
of lfmc~^ has been found for Au + Au collisions at RHIC energies [1] and can be 
confirmed also for Pb + Pb collisions at LHC energies. 

Originally limited to a purely gluonic medium {Nj — 0) the model has been 
extended to include light quarks (Nf = 3). These are considered to be massless 
Boltzmann particles as are the gluons. Binary interactions involving light quarks 
and gluons are computed from leading order pQCD cross sections in small angle 
approximation. Radiative and annihilation processes, gg O ggg, are based on the 
Gunion-Bertsch matrix element [2] 
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2 -H- 3 processes involving light quarks are also computed from (1) based on a 
factorization of the Gunion-Bertsch matrix element into a collisional part and a 
radiation probability I AIgb I = |A^ colli P^- Applying the small angle approximation 
to |Alcoii| , the computation of 2 o 3 processes involving light quarks is reduced to a 
scaling of (1) by color and symmetry factors. 

The cross sections and matrix elements are screened by a Debye mass 
rrijj = dai^cxs J t^^jt?- i^cfg + ^flq) that is dynamically computed from the current 
distributions of gluons and quarks. The Landau-Pomeranchuk-Migdal (LPM) effect 
is modeled via the introduction of a cutoff (A — t) into (1) that effectively discards 
coherent contributions from multiple induced gluon radiation. This is done by a 
comparison of the mean free path A to the formation time of the radiated gluon r. 
More details on the modeling of the LPM effect and on the consequences arising from 
this implementation can be found in [3]. 

3. Nuclear modification factor and elliptic flow for gluons and light quarks 

As established in [4, 5, 6] the matter in BAMPS simulations of Au + Au collisions at 
RHIC energies exhibits a strong degree of collectivity with an integrated V2 that is 
in good agreement with experimental results over a large centrality range for a fixed 
strong coupling of ag — 0.3 and a freeze-out energy density £c = 0.6 GeV fm^'^. These 
parameters are used for all calculations that are presented in this section. 

Using this setup and mini-jet initial conditions withpo = 1.4 GeV, figure la shows 
the nuclear modification factor Raa obtained from BAMPS simulations of central, 
%-10 %, Au + Au collisions at 200 A GeV. The results are both shown on the partonic 
level for gluons and light quarks and on the hadronic level for neutral pions based on 
AKK fragmentation functions [9] . The suppression of high-pT particles in simulations 
with BAMPS is distinctly stronger than the experimentally observed suppression. 
This is due to a) a strong energy loss that is caused by a complex interplay of the 
Gunion-Bertsch matrix element and the effective implementation of the LPM effect 
[3], h) a conversion of quark into gluon jets in 2 — > 3 interactions and c) a small 
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(a) Nuclear modification factor Raa for neutral 
pions from BAMPS simulations of Au + Au at 
200 A GeV compared to PHENIX results [7] at 
%-10 % centrality. Raa of gluons and quarks 
is shown for comparison. Lines indicate Raa 
computed from fits to the parton spectra, while 
symbols indicate Raa computed directly from 
the parton spectra as obtained from BAMPS. 



(b) Nuclear modification factor Raa of charged 
hadrons, gluons and quarks from BAMPS 
simulations of Pb + Pb at fe = fm compared to 
results from ALICE for %-5 % central Pb + Pb 
collisions [8]. Lines and symbols as in figure la. 
For comparison the Raa of charged hadrons 
from simulations of Au + Au at 200 A GeV and 
6 = fm is also shown. 



Figure 1. Nuclear modification factor at RHIC and LHC. 

difference in the energy loss of gluons and quarks caused by the iterative computation 
of interaction rates required by the inclusion of the LPM cutoff into (1). 

Going to the LHC energy of y'SNN = 2.76 TeV and using Pythia initial 
conditions [10], figure lb shows the nuclear modification factor of charged hadrons 
for central Pb + Pb collisions. The results are not significantly different from those at 
RHIC energy and again the suppression of high-py particles is overestimated compared 
to the experimental data. Also the upward trend towards large px present in the 
ALICE data is not reproduced. 

The simulated differential V2 of central and semi-central Pb + Pb collisions at 
2.76 A TeV shows no significant deviation from the Au + Au results at 200 A GeV in 
the low to intermediate px region which is in agreement with experimental findings 
[8]. Following from an increase in the mean transverse momenta, the integrated V2 
however is larger. Figure 2a shows the integrated V2 for Pb + Pb at LHC as a function 
of centrality compared to experimental data from ALICE [8]. Up to roughly 40% 
centrality the agreement is very good, going to more peripheral collisions the simulated 
V2 exhibits a drop that is distinctly more pronounced than in the experimental data. 
The cause of this rapid drop is currently under systematic investigation, comparing 
for example different initial conditions and freeze-out prescriptions. 



4. Heavy quarks 

The BAMPS framework is also extensively applied to the investigation of heavy quark 
phenomena [10, 11, 12]. Figure 2b compares the Raa of Z3-mesons obtained from 
BAMPS to preliminary ALICE data [13]. Heavy quark elastic interactions with 
the gluonic medium {Nj = 0) are implemented with running coupling and a Debye 
screening motivated from HTL calculations [12]. To account for radiative contributions 
the cross section is multiplied by a factor K — A which gives a very good agreement 
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with experimental data for Raa and V2 of heavy flavor electrons at RHIC [12]. It 
will be checked in a forthcoming study whether the implementation of radiative heavy 
quark processes has indeed the same effect as scaling the elastic cross section with a 
constant K factor. 
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(a) Integrated V2 of gluons and light quarks 
{\y\ < 0.8) as a function of centraiity for 
Pb + Pb at 2.76 A TeV compared to the 
measured V2 of charged particles from ALICE 
[14]. 
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(b) Nuclear modification factor Raa for D- 
mesons (|j/| < 0.5) for Pb + Pb collisions at 
2.76 A TeV and b = 4.5 fm for the same setup as 
in [12] . For comparison preliminary data from 
ALICE is shown for 0%-20% centraiity [13]. 



Figure 2. Integrated V2 and D-meson Raa ^t LHC. 



5. Summary 

The transport model BAMPS allows for the simultaneous investigation of both the 
suppression of high-pT particles and the collectivity of bulk particles consistently 
within a common setup that handles gluons as well as light and heavy quarks. 
Employing the same parameters that give a good description of the integrated elliptic 
flow at RHIC energies, the suppression of high-p^ particles from gluons and light 
quarks is too strong. The same holds for simulations of Pb + Pb collisions at LHC, 
where no significant differences in the suppression pattern and the differential elliptic 
flow are observed compared to RHIC simulations. The integrated elliptic flow at LHC 
is well described up to 40 % centraiity, for more peripheral collisions the simulated V2 
drops too fast. Fixing the interaction of heavy quarks via elastic interactions with 
the medium to the suppression of heavy flavor electrons at RHIC by means oi a K 
factor, the elliptic flow of heavy flavor electrons at RHIC as well as the suppression of 
D-mesons at LHC can be reproduced. 

Future studies will investigate the important effect of a running coupling on the 
Raa of gluons and light quarks and also systematically explore the implementation 
of the LPIvI effect. These modifications are qualitatively expected to bring the results 
for the nuclear modification factors into better agreement with experimental data. 
The consequences of radiative processes on the dynamics of heavy quarks will also be 
investigated in an upcoming study. 
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